Measuring Flavor Ratios of High-Energy Astrophysical Neutrinos 



John F. Beacom/E Nicole F. Bell,i^ ^-Q Dan Hooper,3^[3 Sandip Pakvasa,^'^'! and Thomas J. Weiler^-^'^ 

' NASA/Fermilab Astrophysics Center, Fermi National Accelerator Laboratory, Batavia, Illinois 60510-0500 
'^Kavli Institute for Theoretical Physics, University of California, Santa Barbara, California 93106 
'^Department of Physics, University of Wisconsin, Madison, Wisconsin 53706 
^Department of Physics and Astronomy, University of Hawaii, Honolulu, Hawaii 96822 
"^Department of Physics and Astronomy, Vanderbilt University, Nashville, Tennessee 37235 

(Dated: July 1, 2003) 

We discuss the prospects for next generation neutrino telescopes, such as IceCube, to measure the 
flavor ratios of high-energy astrophysical neutrinos. The expected flavor ratios at the sources are 
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current, and all flavors neutral-current), muon tracks {v^ charged-current only), and tau lepton 
lollipops and double-bangs {vr charged-current only). The peak sensitivities for these interactions 
are at different neutrino energies, but the flavor ratios can be reliably connected by a reasonable 
measurement of the spectrum shape. Measurement of the astrophysical neutrino flavor ratios tests 
the assumed production mechanism and also provides a very long baseline test of a number of exotic 
scenarios, including neutrino decay, CPT violation, and small-^m"^ oscillations to sterile neutrinos. 
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I. INTRODUCTION 

A new generation of detectors proposed or already un- 
der construction will have the sensitivity to open a new 
window on the universe in the form of high-energy neu- 
trino astronomy. Besides probing distant and mysterious 
astrophysical sources, these data will offer unprecedented 
sensitivity for testing fundamental neutrino properties. 

Under-ice and under-water optical Cerenkov detec- 
tors are sensitive at roughly TeV-PeV neutrino ener- 
gies. The AMANDA detector Q at the South Pole 
has already been taking data for several years, as has 
the smaller Lake Baikal detector 0. First results from 
the AMANDA-BIO phase have been reported [J 
and results from several years of running in the larger 
AMANDA-II configuration 1| are eagerly awaited. Con- 
struction of the much larger (km'^) IceCube detector 5] at 
the South Pole is scheduled to begin later this year. The 
ANTARES and NESTOR 01 collaborations, draw- 
ing on the lessons learned by the DUMAND project 
are currently deploying their detectors in the Mediter- 
ranean Sea. While smaller than IceCube, these detectors 
will have lower thresholds due to denser instrumentation. 
Another Mediterranean detector, NEMO has been 
proposed. These detectors operate under water or ice to 
shield the atmospheric muon flux; since they primarily 
look for upgoing neutrino-induced events, the detectors 
in the southern and northern hemispheres offer comple- 
mentary views of the sky. 
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Detectors making use of the Earth's atmosphere as a 
target volume are sensitive at roughly EeV-ZeV neutrino 
energies. Fly's Eye [13 and AGASA designed to 
detect atmospheric showers induced by cosmic ray pro- 
tons, are also sensitive to ultra-high energy neutrino pri- 
maries if they consider penetrating horizontal showers, 
i.e., events that initiate only at great slant depth in the 
atmosphere. The bounds on ultra-high energy neutrino 
fluxes are reviewed in Ref. [l^]- The HiRes detector (an 
upgrade of Fly's Eye) is running '13?|, and the Pierre 
Auger detector scheduled to be completed in 2005, 
is already taking data with their engineering array. Both 
should significantly improve the sensitivity to ultra-high 
energy astrophysical neutrinos. Orbiting detectors will 
soon provide a new window on neutrinos as well. The 
EUSO experiment 0| is scheduled for deployment on 
the International Space Station in early 2009, and it may 
evolve into a larger satellite mission named OWL [l6j . 

Radio Cerenkov detectors are sensitive at roughly PeV- 
ZeV neutrino energies. As pointed out long ago by 
Askaryan a Cerenkov signal is proportional to the 
square of the net charge of the shower within a wave- 
length. In turn, the net charge is roughly proportional to 
the shower energy. Thus, the long- wavelength radio sig- 
nal rises as the neutrino energy squared, instead of just 
the neutrino energy, as for other techniques. The RICE 
antennas, co-deployed with AMANDA, have recorded 
first results Recently, NASA approved a very in- 

teresting experiment, ANITA |T3, which will carry radio 
antennas on a balloon above the South Pole in an attempt 
to detect the long-wavelength tail of sub-ice showers from 
Earth-skimming neutrinos [23|. The GLUE experiment 
has reported limits based on the non-observation of a ra- 
dio signal from neutrino interactions in the surface of the 
Moon ||2l|. The proposed SALSA detector could use 
large salt domes as a radio Cerenkov medium. 
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Although the goals of these experiments are common, 
the detection strategies and systematic issues are not. 
Each will infer the direction of arriving events, thereby 
enabling point-source neutrino astronomy. Each will 
measure the neutrino flux over some energy range. One 
purpose of this paper is to ask how well the primary neu- 
trino spectrum can be reconstructed from the observed 
spectral information. The primary spectrum reveals the 
dynamics of the cosmic engine. Reviews of high-energy 
neutrino astronomy are listed in Ref. |23j |. 

Besides energy and direction, an additional piece of in- 
formation is carried by arriving neutrinos: flavor. Neutri- 
nos are known to come in three flavors, electron-, muon-, 
and tau-type. The bulk of this paper focuses on what we 
can learn from flavor identification of the incoming neu- 
trinos, and whether flavor-tagging is feasible with pro- 
posed experiments. We consider the IceCube experiment 
in detail due to its large effective volume and ability to 
observe showers, muon tracks and events unique to tau 
neutrinos. We emphasize that other experiments will also 
have some capabilities to discriminate among neutrino 
flavors. 

In the next section, we look at theoretical motivations 
for flavor discrimination. In the sections after, we look 
at inferences of the incident neutrino spectrum, and the 
neutrino flavors from various event signatures. 



II. WHY NEUTRINO FLAVOR 
IDENTIFICATION IS INTERESTING 

Neutrinos from astrophysical sources are expected to 
arise dominantly from the decays of charged pions (and 
kaons) and their muon daughters, which results in ini- 
tial flavor ratios, (piy^ : 4)^^ : (j),^^, of nearly 1:2:0. 
The fluxes of each mass eigenstate are then given by 
(jjj = Y.a (t>lT'"'\Uaj\^ ■ The Ucj are elements of the 
neutrino mass-to-flavor mixing matrix, defined by Iva) = 
"^jUajli^j)- The propagating mass eigenstates acquire 
relative phases, giving rise to flavor oscillations. How- 
ever, these relative phases are lost, since 6m^ x L/E ^ 1, 
and hence uncertainties in the distance L and the energy 
E will wash out the relative phases. Thus the neutri- 
nos arriving at Earth are an incoherent mixture of mass 
eigenstates with the proportions given above. 

For three neutrino species, as we assume throughout, 
there is now strong evidence from atmospheric and re- 
actor neutrino data suggesting that and Vr are max- 
imally mixed and Ues is nearly zero. This twin happen- 
stance (^atm — 45° and Ues = 0) leads to two remarkable 
conclusions. The first is that each mass eigenstates con- 
tains an equal fraction of and Vt- The second is that 
in the mass eigenstate basis, the neutrinos are produced 
in the ratios 1:1:1, independent of the solar mixing 
angle, and thus arrive at Earth as an incoherent mixture 
of mass eigenstates with these same ratios. This implies 
democracy in the detected flavor ratios as well, since since 
= 1 in any basis. 



So there is a fairly robust prediction of 1 : 1 : 1 flavor 
ratios for measurements of astrophysical neutrinos 
I25I I . The first task of flavor measurement is to check this 
prediction. Could the flavor ratios differ from 1:1:1? 
The astrophysics could be different than is outlined here. 
For example, if the charged pion decays promptly but 
the daughter muon loses energy before decaying, then at 
high energy the flux may approximate 4>^^ : (j),^^ : (f>^^ 
0:1:0 [23|. This leads to mass eigenstate ratios of 
^1:2:3, and measured flavor ratios of ~ 1 : 2 : 2. 

Alternatively, the particle physics could be different 
than assumed here (three stable neutrinos). For exam- 
ple, the heavier mass eigenstates could decay en route to 
Earth. This leads to markedly different detected flavor 
ratios, as extreme as either 6 : 1 : 1 or : 1 : 1 for the 
normal and inverted hierarchies, respectively j27l |. Also, 
Barenboim and Quigg have pointed out that CPT viola- 
tion in neutrino mixing could also lead to anomalous fla- 
vor ratios 28j. Finally, neutrinos could be pseudo-Dirac 
states, in which case the three active neutrinos have ster- 
ile partners, with which they are maximally mixed with 
very tiny Sm^ splittings; these oscillations might only be 
effective over cosmological distances ^29] . 

In all cases investigated, the — Vt symmetry ensures 
that and i/,- arrive at Earth in equal numbers. Given 
how robust that prediction is, and that Vt is the most 
difficult flavor to identify in IceCube, we will mostly fo- 
cus on the (l>u^ : 0^ ratio. We will show that this can be 
determined by measuring the rates of shower and track 
events in IceCube. Additionally, we discuss the identifi- 
cation of Vr ; though the expected yields are very small, 
the signals are very distinctive, and the detection of even 
a single event would confirm the presence of the flux. 



III. BASICS FIRST: MEASURING THE 
NEUTRINO SPECTRUM WITH MUONS 

For under-ice or under-water detectors, muon events 
provide the most useful signal from which to infer the 
neutrino spectrum. This is especially true at lower en- 
ergies (100 GeV-several TeV) due to the higher energy 
threshold for showers. In this section we explore recon- 
struction of the neutrino spectrum from observed muon 
events in an IceCube-type detector. There are challenges 
in doing this, however. 

First, muons can be created in interactions far from the 
detector, and lose a considerable fraction of their energy 
before being measured. This problem can be circum- 
vented, however, by considering only muon events with a 
contained vertex, in which the muon track begins within 
the detector volume. At energies near or below the TeV 
scale, many of the observed muons will have contained 
vertices. At higher energies, when the range of muons 
is considerably longer, fewer of the resulting events will 
have this feature. 

Second, separating the atmospheric neutrino back- 
ground from any astrophysical neutrino signal can be 
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FIG. 1: The distribution of observed muon energies for a 
neutrino spectrum of El^^dN^^/dE^^ — 10~^ GeV cm~^ s~^ 
for one year. The dashed hue is for a flux at a horizontal 
zenith angle, and the dotted line for an upgoing flux (through 
the Earth). Note that the muon energy at production may be 
considerably larger than that observed due to the long muon 
range. The error bars are slightly ofl'set for clarity. 



FIG. 2: The distribution of observed contained vertex 
muon energies for a neutrino spectrum of El^^dN^^J dE^^^ = 
10-'' GeV cm s for one year. The dashed line is for a 
flux at a horizontal zenith angle, and the dotted line for an 
upgoing flux (through the Earth). The error bars are slightly 
offset for clarity. 



difRcult. The atmospheric neutrino spectrum is well- 
modclcd, and so in principle can be subtracted from the 
data. At energies at or above about 100 TeV, the as- 
trophysical neutrino flux is likely to be above the more 
steeply falling atmospheric neutrino flux. At more mod- 
est energies, the angular and temporal resolution of a 
neutrino telescope will be needed to effectively remove 
backgrounds. The atmospheric neutrino flux at 1 TeV 
(dominated by + is ~ 10"^ GeV cm'^ s"^ in 
a (1 deg)^ bin, where the bin size was chosen to reflect 
the angular resolution for these events. Especially for the 
lower flux we consider below, identification of astrophys- 
ical sources may require temporal information as well. 
For example, gamma-ray bursts typically have durations 
on the order of seconds. Taking the known catalogs of 
gamma-ray bursts as a guide, about 10 events per square 
kilometer per year are expected. Using timing and di- 
rectional information, these events are essentially back- 
ground free [sij . Tau neutrino detection for gamma-ray 
bursts has been studied in Ref . j32l | . A similar technique 
could be used for blazars or other transient sources. 



A. Muon Tracks 

After a high energy muon is produced, it undergoes 
continuous energy loss as it propagates, given by 



dE 
dX 



where a = 2.0 MeV cm^/g and /3 = 4.2 x 10"*^ cm^/g [32 
The muon range is then 



a + l3En 



thr 



(3.2) 



-a- PE , 



(3.1) 



where -E*'" is the minimum muon energy triggering the 
detector. Typically, E^^' - 50 - 100 GeV for deep 
ice or water detectors. Above an energy of a TeV, the 
muon range rises logarithmically as ^ lii{E / E^^" ) times 
2.4 km water equivalent. Since this typically exceeds the 
size of the detector, the muon energy cannot be mea- 
sured by the muon range. Further, since the muons are 
always fully relativistic, the Cerenkov angle and intensity 
are constant and thus cannot be used to infer the muon 
energy. However, the muon energy inside the detector 
can be inferred by the rate of energy deposition in the 
form of showers from catastrophic bremsstrahlung [5ll33|. 
The muon range is substantially less than the muon decay 
length; for an illustration of the length scales for neutrino 
interactions, mu and tau range, and mu and tau decay, 
see Fig. 1 of Ref. [13 . 

The energy of the muon faithfully represents the neu- 
trino energy since the charged-current differential cross 
section is strongly peaked at y = \ — E^/Ei, ~ 0, and 
{y) ~ 0.2 [s^. The kinematical angle of the muon rel- 
ative to the neutrino direction is about 1°/ ^jE„/\ TeV, 
and the reconstruction error on the muon direction is on 
the order of 1°. 

The probability of detecting a muon neutrino traveling 
through the detector via a charged-current interaction is 
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FIG. 3: The distribution of observed muon energies for a 
neutrino spectrum of El^^dN^^/dE^^ — 10~^ GeV cm~^ s~^ 
for one year (dashed), as in Figs. and |21 compared against 
spectra proportional to E^^'"^ (solid) and E^^'^ (dotted), nor- 
malized to the same number of events in the first energy bin. 
All rates are for a horizontal zenith angle source. Note that 
the muon energy at production may be considerably larger 
than that observed. The error bars are slightly offset for clar- 
ity. 

then given by 

P,^^^ ~ pNAaRf, , (3.3) 

where p is the target nucleon density, Na is Avogadro's 
number, and a is the neutrino-nucleon total cross sec- 
tion 36J. 

B. Spectral Results 

In Figs, n and [21 we show the distribution of ob- 
served muon energies for a muon neutrino spectrum of 
El^dN^JdE^^ = 10^"^ GeV cm-^ s"! for one year (or 
El dN^JdE^,^ = 10-8 GeVcm-2 for ten years). 
The dashed Hues correspond to the flux at a horizon- 
tal zenith angle, and the dotted lines to a upgoing flux 
(through the Earth). Detection prospects for horizon- 
tal neutrinos are enhanced from long distances of ice in 
which muons can be produced. Upgoing neutrinos, al- 
though also with this advantage, can be absorbed in the 
Earth, degrading their event rate. A distribution of neu- 
trino sources over the sky produces a spectrum in be- 
tween these two extreme curves. Fig. ^ shows the spec- 
trum of all observed muons, while Fig.|2shows only muon 
events with a contained vertex. The events are binned by 
energy decade, with 68% confidence levels shown. The 
dashed and dotted lines, shown for comparison, are the 
results of our Monte Carlo with much greater statistics. 



FIG. 4: The distribution of observed contained vertex 
muon energies for a neutrino spectrum of E^^dN^^J dE^^ = 
IQ-^ GeV cm"^ s"^ for one year (dashed), as in Figs. 0and 
121 compared against spectra proportional to E~'^'^ (solid) and 
E~^'^ (dotted), normalized to the same number of events in 
the first energy bin. All rates are for a horizontal zenith angle 
source. The error bars are slightly offset for clarity. 



Note that the sizes of the energy bins was selected for 
statistical purposes. The energy resolution of neutrino 
telescopes is considerably more precise at these energies. 

In Figs. and 21 we compare the observable horizon- 
tal muon spectra resulting from neutrino spectra propor- 
tional to £'-2 (dashed), E-^-^ (solid), and E-'^-^ (dot- 
ted), normalized to the same number of events in the 
first energy bin. By comparing the numbers of events 
in adjacent bins. Fig. O demonstrates that for a single 
power-law flux (and our choice of normalization), the 
spectral slope can be determined to approximately ten 
percent up to tens or perhaps hundreds of TeV. Fig. ^ 
demonstrates the ability to make such as measurement 
with only muons with contained vertices. Comparing 
Figs. and ^ makes it clear that limiting the data to 
muons with contained vertices alone weakens the ability 
to resolve similar spectral slopes, especially at higher en- 
ergies. This is because of the statistical limitations. How- 
ever, the fact that muons with contained vertices yield 
the total muon energy, whereas throughgoing muons of- 
fer only lower limits to the energy, has its virtue. The 
contained-muons measurement is more useful for the case 
of a strongly broken power-law, such as is predicted in 
the case of gamma-ray bursts. Typically, the gamma-ray 
burst neutrino spectrum is expected to follow a E^^ spec- 
trum up to some break energy on the order of hundreds 
of TeV. Above this energy, it steepens to a spectrum 
proportional to E~^. For such a spectrum, the break 
could be observed with the contained vertex muons mea- 
surement, and the slopes could be then measured more 
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precisely using all observed muons. 

We note that the flux we employ for illustration, 
El^dN^JdE^^ = 10-^ GeV cm-2 s-\ is on the or- 
der of the Waxman-Bahcall bound [s^- The Waxman- 
Bahcall bound pertains to the diffuse neutrino flux from 
sources optically thin to protons, normalized to the mea- 
sured cosmic-ray flux at ^ eV. As emphasized by 
Mannheim, Protheroe and Rachen [s^ , it does not apply 
for optically thick, or "hidden" sources [s^j for galactic 
sources, or for sources not emitting protons at energies 
~ 10^^ eV. Microquasars and supernovae remnants pro- 
vide examples of both of the latter categories. The cur- 
rent experimental limit on the high-energy neutrino flux, 
from the AMANDA experiment, is more than an order of 
magnitude larger than the Waxman-Bahcall bound 
IceCube is expected to reach well below the Waxman- 
Bahcall bound . 

Nature's flux could be larger than we assume here. It 
could also be smaller. If it is smaller, then integration 
times larger than the one year we assume here are needed 
to compensate. Furthermore, larger detectors, such as 
an extension of IceCube, are likely to be constructed in 
the future, making more conservative choices of the neu- 
trino flux easier to study. To reduce backgrounds, only 
those events associated with known sources should be 
considered. In that case, the remaining flux of neutri- 
nos which could be used for such a study may be sig- 
nificantly reduced. Given this consideration, perhaps a 
choice of dN^/dE^,, = 10"® GeV cm'^ s"^ over ten 
years, which we also discuss, could be considered a more 
realistic choice. 



IV. FLAVOR IDENTIFICATION 

Although the ratios of neutrino flavors are not directly 
measurable, they can be inferred from complementary 
classes of events in neutrino telescopes. In this sec- 
tion, we restrict our attention to IceCube, which will be 
capable of identifying showers from both charged- and 
neutral-current events, muon tracks, and certain tau neu- 
trino events. Muon tracks have been discussed above. 
The probabilities for detecting the different neutrino fla- 
vors are illustrated in Fig. |5l and will be discussed in 
detail now. 

Electron and muon neutrinos above about 100 TeV 
are absorbed in Earth by their charged-current interac- 
tions. Tau neutrinos also interact, but regenerate i^r by 
the prompt decays of tau leptons these decays also 
produce a secondary flux of Vf, and |35l l4lj| . 



A. Showers 

All neutrino flavors undergo an identical neutral- 
current interaction producing a hadronic shower and 
nothing else. The shower energy typically underestimates 




K (GeV) 

FIG. 5: Probabilities of detecting different flavors of neutri- 
nos in IceCube versus neutrino energy, described in detail in 
the text. The upper solid line is the probability of a hori- 
zontal creating a detectable muon track, and the dashed 
line is for downgoing Vfj,. The dotted line is the probability 
for ly^ to create a detectable shower (above 1 TeV), consid- 
ering both charged-current and neutral-current interactions; 
the kink occurs when the neutral-current showers come above 
threshold. The dot-dashed lines are the probabilities for Vt to 
make lollipop events (upper) and double-bang events (lower). 



the neutrino energy by a factor ranging from ~ 3 around 
a TeV, to a factor of ~ 4 at an EeV [Sg ■ 

A charged-current interaction of produces an elec- 
tron that immediately creates an electromagnetic shower. 
If the electromagnetic shower is measured with the 
hadronic one, the total shower energy is the incident 
Ve energy. In principle, electromagnetic and hadronic 
showers are distinguishable by their respective muon con- 
tent, absent for electromagnetic and present for hadronic 
showers. We will not assume that these can be distin- 
guished, as it is expected to be very difficult. 

In a charged-current interaction, the muon track al- 
ways emerges from the shower, because of the long muon 
range, so these events don't contribute to the shower rate. 

The Vr charged-current interaction produces a 
hadronic shower and a tau track. The tau decay path- 
length is 7cr ~ 50 {Er/ PeV) m. Below a few PeV, the 
tau track is too short to be separated from the shower, 
and so these events will contribute to the shower rate. 
At higher energies, the tau track will extend beyond the 
initial shower, and then the tau will decay to produce a 
second shower. This creates identifiable double-bang and 
lollipop events, discussed below. 

Showers are seen by the detector as photoelectrons dis- 
tributed over a ~ 100 m radius sphere for a TeV shower 
(~ 300 m radius for showers with PeV energies). The 
shower must at least be partially contained within the de- 
tector volume in order to be detected. Since shower sizes 
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are relatively small compared to muon ranges, the effec- 
tive volume for these events is substantially less than for 
charged-current f^, interactions. Also, the energy thresh- 
old for showers is generally larger than for muon tracks. 

The probability of detecting a neutrino by a shower 
produced by a neutral-current interaction is given by 



V — ^shower 



pNaL 



da 

-i-dy 



(4.1) 



where a is the neutrino- nucleon cross section j3fil |. y is 
the energy fraction transferred from the initial neutrino 
to the hadronic shower, and L is the length of the detec- 
tor. For charged-current electron neutrino interactions, 
however, the additional electromagnetic shower means 
that all of the neutrino energy goes into the shower, and 



P, 



'y^showor — pNaCtL . 



(4.2) 



A similar treatment is used for charged-current tau neu- 
trino interactions which do not produce a double-bang 
or lollipop event (see below). In real experiments, the 
shower threshold is not a step function as we adopt here. 
IceCube will have shower energy resolution of about ±0.1 
on a logj^Q scale and will be able to reconstruct the neu- 
trino direction to about 25°. 



B. Double-bang and lollipop events 

Double-bang and lollipop events are signatures unique 
to tau neutrinos, made possible by the fact that tau lep- 
tons decay before they lose a significant fraction of their 
energy [s^]. Double-bang events O 43] consist of 
a hadronic shower initiated by a charged-current inter- 
action of the Pr followed by a second energetic shower 
(hadronic or electromagnetic) from the decay of the re- 
sulting tau lepton. Lollipop events consist of the second 
of the two double-bang showers along with the recon- 
structed tau lepton track (the first bang may be detected 
or not). Inverted lollipops, consisting of the first of the 
two double-bang showers along with the tau lepton track, 
are not as useful as they will often be confused with a 
hadronic shower in which a ~ 100 GeV muon is pro- 
duced (because of the higher lepton mass, tau tracks suf- 
fer much less catastrophic bremsstrahlung than muons at 
the same energy). We do not consider inverted lollipops 
for this reason. 

The range of a tau lepton is bounded from above by 
its lifetime in the lab frame. A tau lepton of energy Er 
has a mean lifetime given by 

RT{Ey^,y)^ CTr^ CTr, (4.3) 

mr rrir 

where m,- and Tt are the mass and rest-frame lifetime. 

Following Ref. [23. li^ . |43| , the conditions which must 
be fulfilled for the detection of a double-bang event are: 



• The tau neutrino must interact via the charged- 
current, producing a hadronic shower of sufficient 
energy to trigger the detector (~1 TeV) inside of 
the detector volume. 

• The tau lepton produced in the interaction must de- 
cay inside the detector volume, producing an elec- 
tromagnetic or hadronic shower of sufficient energy 
to trigger the detector. 

• The tau lepton must travel far enough such that the 
two showers are sufficiently separated to be distin- 
guished from each other. 

At the energies required for the third condition to be 
satisfied, both the showers will be energetic enough to 
easily fuUfil the threshold requirements. The probability 
for a double-bang event, per incident tau neutrino, is 



Pdb(S.J 



pNa / dy— 
Jo dy 

pNa / dy^ {L 
Jo dy 



dx 



Rr 



Rr)e 



RrC 



(4.4) 



The track integration from Xmin to L includes tau lengths 
that are larger than Xmin (for shower separation) and 
smaller than L (so that both showers are contained in 
the detector.) The exponential in the integral samples 
over the decay length distribution. 

Conditions which must be fulfilled for the detection of 
a lollipop event are: 

• The shower produced by the decay of the tau lepton 
must occur within the detector volume and be of 
sufficient energy to trigger the detector (~ 1 TeV) . 
This may be a hadronic or electromagnetic shower. 

• The track of the tau lepton must be long enough 
(within the detector) to be reconstructed and sep- 
arable from the shower. For photomultiplier tube 
spacing of ~ 125 meters (IceCube's horizontal spac- 
ing is 125 meters, their vertical spacing is a signifi- 
cantly smaller 17 meters), reasonable values for the 
minimum tau range, Xmim are 200 — 400 meters. 

The probability for a lollipop event, per incident tau neu- 
trino, is 



f da 

PioUipop{E^^) ^ pNa{L - Xmin) / dy—e' 

Jo dy 



(4.5) 

Note that the energy threshold for both double-bangs 
and lollipops, resulting from the requirement Rr ^ a;i„in, 
IS given by E^J^"^ ~ 5 (xmin/250 m) PeV. For double-bang 
events, the energy threshold for the first shower, Egh ~ 
y Eu^ ^ TeV puts a lower limit on the y-integration, but 
negligibly so for PeV scale neutrino energies. 
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bangs separated by distances larger up to a factor -s/S- 
The analytic expressions and the Monte Carlo calculation 
shown in Fig. El are in good agreement, a considerable im- 
provement over the approximations given in Rcf . "4^, '43| . 
In particular, the important threshold region is charac- 
terized very well. The remaining minor discrepancies we 
attribute to the low statistics of the Monte Carlo; our 
formulas should be quite accurate. In addition, the ap- 
proximation in Eq. (|4.6|l is seen to yield quite accurate 
results (which indicates that the variation of cross section 
with y is relatively unimportant in these calculations). 
Thus we have good confidence that Eqs. (|4.4I) to l|4.8|) 
may join with Eqs. 1)3. 3() to (|4.2() to complete the set of 
probabilities of event topologies. 



FIG. 6: Probability of observing a lollipop (upper) or double- 
bang (lower) event in IceCube per incident tan neutrino. The 
points (X's and O's) represent Monte Carlo results while the 
solid lines represents the analytic expressions of Eqs. (14.411 and 
14.511 and the dashed lines represents the approximations of 
Eqs. I|4.7fl and 114. 8|l . We have set L = 1 km and a;min = 250 
m. Above about 20 PeV, the tau track length becomes long 
enough that it may exceed the thickness of the ice for near- 
downgoing events, introducing zenith-angle dependence. 



The expressions in Eqs. 14.4|l and H4.5|l can be further 
simplified with the approximation 



da 
dy 



(4.6) 



where {y) ~ 0.25 at PeV scale energies and a is the total 
cross section. Thus, 



Pdb(£;.J ^ pNa<j 



flonipop(£'^,) — PN ACT {L - Xmin) 



Rr) e-^ 

v=(y) 



d„/-Rx 



(4.7) 



(4.8) 

y={y} 



Fig. El shows the lollipop and double-bang detection 
probabilities. Note that near threshold, the majority of 
lollipop events are also double-bang events. At higher 
energies, the size of the detector excludes many double- 
bang events but not lollipops. 

The analytic expressions in Eqs. H4.4|l to H4.8|l are orig- 
inal to this work. To check their validity, we constructed 
a Monte Carlo to calculate the probabilities of observing 
double-bang or lollipop events in IceCube. This simu- 
lation took into account the distribution of y values in 
charged-current interactions as well as the variation in 
decay lengths. In common with the analytic results, the 
Monte Carlo treats the detector one dimensionally. We 
note that this is a conservative approximation as taus 
that travel diagonally across the detector could have 



C. Determining the Muon/Shower Ratio 



The spectrum of shower events is more difficult to infer 
from data than is the muon spectrum due to lower statis- 
tics, but it may be possible (i^. Relative disadvantages 
are the considerably higher shower energy threshold [3 
and the fact that muon event rates benefit from long 
muon ranges. However, in contrast to the muon signals, 
shower events are produced by all flavors of neutrinos, 
and shower energies more faithfully represent the neu- 
trino energy than do muon tracks. We do not attempt to 
relate the observed shower spectrum to the spectrum of 
incident neutrinos. Rather, we assume that, due to os- 
cillations, the neutrino spectrum shape is independent of 
flavor; i.e., the neutrino spectrum inferred from a mea- 
surement of the muon spectral shape is universal. We 
may use this universal spectrum to produce shower events 
and muon tracks, and then compare the total number of 
shower events to the total number of muon events to 
obtain the flavor ratios. Recall that i^^ — symmetry 
means 4>^^ : 4>^^ =1:1, and so two independent observ- 
ables are sufficient to determine all three fiavor ratios. 

For a neutrino spectrum of E^^dNi^^/dE^^ — 
10^^ GeV cm^^ s^^ for one year, we expect, on average, 
323 muon events (186 of which have contained vertices). 
Given an equal number of each neutrino flavor, we predict 
only 36 shower events. At the 68% confidence level, this 
allows for a measurement of A^muons /-^showers = Q.Oti'g 
for horizontal sources, about a 20% uncertainty. For 
upgoing sources, a measurement of A^muons /^showers = 
8.5ii[g results similarly. The difference comes from the 
fact that upgoing neutrinos are absorbed at high energies, 
where muon events are more likely to occur, thus slightly 
lowering the muon to shower ratio for upgoing events rel- 
ative to horizontal events. These results use the natural 
energy thresholds of the detector, ~ 100 GeV for muons 
and 1 TeV for showers. 
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FIG. 7: The relationship between the muon-to-shower ra- 
tio and the Ue fraction, assuming j/^ — i/t symmetry and an 
power law spectrum. The muon energy threshold is 100 
GeV, and the shower threshold is 1 TeV. Horizontal sources 
are assumed. The solid line is the central predicted value. 
The dashed lines represent the 68% confidence interval for a 
spectrum of El^dN^^/dE^^ = 10~^ GeV cm s for one 
year. The dotted lines represent the 68% confidence interval 
for a spectrum five times smaller. 



D. Signatures Unique To Tau Neutrinos 

Given the flux we have considered in this paper, 
E^^dNy^/dE^^ = 10~^ GeV cm~^ s^^ for one year, we 
predict on the order of a 50% chance of observing a lol- 
lipop event, and similarly for a double-bang event. Thus 
IceCube is unlikely to provide a stringent probe of the 
tau neutrino flux. The lack of observed double-bangs 
or lollipops in IceCube would not reveal much, though 
the positive identiflcation of even a single such event by 
IceCube would indicate the important existence of a tau 
neutrino flux on the order of the flux we consider here. 

Even larger detectors are needed to exploit the double- 
bang and lollipop features. With several of these events, 
flavor ratios could be easily reconstructed, and the i/^ — i^^- 
symmetry tested. Since even larger detectors are even 
farther into the future, we do not consider the double- 
bang and lollipop signatures further in this work. 



V. INFERRING NEUTRINO FLAVOR RATIOS 

Although the neutrino flavor ratios are not directly ac- 
cessible at neutrino telescopes, the indirect flavor infor- 
mation collected from such experiments, i.e., the ratios of 
muon, shower and tau-unique events, can be very useful 
in inferring flavor information. 

In the previous section, we showed that for a neutrino 
spectrum of E'^, dN^ IdE^, = 10"'' GeV cni"^ s"^ for 

f_L M ' M 



0.0 0.2 0.4 0.6 0.8 1. 

Ji'raction 

FIG. 8: Same as Fig. |7| except that a 1 TeV muon energy 
threshold has been imposed to reduce the effect of the uncer- 
tainty in the muon spectrum. 



one year, IceCube could determine the ratio of the muon 
events to shower events with uncertainties on the order of 
20%. Given — Ur symmetry, this ratio can be used to 
deduce the ratio of electron neutrinos to either muon or 
tau neutrinos. Clearly, precision measurements of these 
quantities are unlikely to be determined in this fashion. 
Fortunately, in some interesting theoretical scenarios, the 
predictions for deviations from a 1:1:1 flavor ratio are so 
extreme as to not require greater precision. The neutrino 
decay model |27j provides a splendid example of possibly 
large flavor deviations. 

In Fig. Q we show the relationship between the muon- 
to-shower ratio and Ue fraction, assuming i^^ — Vr sym- 
metry and an E^"^ power law spectrum. We define the 
fraction as the fraction of neutrinos with electron flavor 
(1 for all electron neutrinos, for no electron neutrinos). 
The solid line is the central predicted value. The dashed 
lines represent the 68% confidence interval for our chosen 
spectrum. The dotted lines represent the 68% confidence 
interval for an exposure five times smaller. Note that the 
Ve fraction of 1/3, expected from known particle physics, 
can be measured to a range of 0.26-0.37 or 0.18-0.44 for 
the two fluxes shown in Fig. |7| respectively. These mea- 
surements are sufficient to test interesting speculations 
such as neutrino decay. 

If the neutrino spectrum is not well measured, then 
the relationship between the flavor ratios and the 
muon/shower ratio will be altered. In the case of poor 
muon spectrum resolution, some steps could be taken to 
reduce the related uncertainties. For example, the num- 
ber of shower events could be compared to muon events 
above a threshold near 1 TeV, rather than the experi- 
mental threshold near 100 GeV. By doing this, the effect 
of the spectral slope is reduced, as shower events and 
muon events only from neutrinos with similar energies 
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TABLE I: Summary of the muons/showers ratios expected 
for selected scenarios. The decay scenarios with normal and 
inverted neutrino mass hierarchies are taken from Ref. 



ratios at 
source 


decays 


ratios at 
Earth 


fraction 


muons /showers 
(Fig. 7 central value) 


1:2:0 


none 


1:1:1 


0.33 
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Normal 


6:1:1 


0.75 


1.5 




Inverted 


0:1:1 





40 


0:1:0 


none 


1:2:2 


0.2 


14 



are compared. In Fig. |S1 we show results analogous to 
Fig. □ but with a muon energy threshold of 1 TeV im- 
posed (Dutta, Reno, and Sarcevic considered a show- 
ers/muons ratio to test three-flavor active neutrino oscil- 
lations against no oscillations, assuming standard flavor 
ratios at production; our results are in reasonable agree- 
ment at the one point of common consideration, a 
fraction of 1/3). With this choice of threshold, a mea- 
surement ofQ; = 2.0±0.2 corresponds to an uncertainty 
on the order of only 20% in the predicted muon/shower 
ratio. A v^, fraction of 1/3 can be measured to a range of 
0.22-0.42 or 0.09-0.61 for the two fluxes shown in Fig.il 
without using any information from muons below 1 TeV. 
With Fig. [7| as our guide, we summarize in Table d the 
muon-to-shower ratio expected for each of various astro- 
physical neutrino models. 



VI. CONCLUSIONS 

Next generation high energy neutrino telescopes, such 
as the kilometer scale experiment IceCube, will be ca- 
pable of observing muon tracks produced by muon neu- 
trinos, shower events from all flavors of neutrinos, and 
possibly the double-bang and lollipop topologies unique 
to tau neutrinos. Using these features, it will be possible 



to infer the flavor ratios of astrophysical neutrino fluxes. 
These can be used to probe properties of the sources 
and to test for new physics beyond the standard model 
of neutrino physics |27L 123 . l29l | . 

Muon events, being the most numerous in under- 
ice/ water detectors, are the most useful for measuring the 
astrophysical neutrino spectrum. We demonstrate that 
this can be accomplished with acceptable accuracy, given 
a sufflciently large but realistic neutrino flux. Knowledge 
of the spectral shape then allows us to make comparisons 
of the total number of event types (muons, showers, or 
tau-unique events) to infer the neutrino flavor ratios. 

Assuming t'^ — Vr symmetry as indicated by oscillation 
data, the ratio of showers to muons provides sufficient 
information to determine all three neutrino flavor ratios. 
Large deviations from the flavor ratios predicted by os- 
cillations (1:1:1) will likely be observable. Furthermore, 
tau-unique events (double-bangs and lollipops) provide 
another tool with which to address flavor identification; 
however, unless there is a larger than expected tau neu- 
trino flux, such events will probably be too small in num- 
ber in kilometer-scale detectors to provide useful informa- 
tion beyond the existence of astrophysical tau neutrinos. 
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(Erratum Date: 15 June 2005) 

Due to a programming error, some of the results presented in our published paper [1] are incorrect. Here we revise 
Figs. 1-4,7,8 and some of the associated discussion, including the final column of Table I. The conclusions of the paper 
are mostly unchanged, and in some cases improved, despite large changes in some of the intermediate results. 

Figures 1-4 show the observable muon spectra for the IceCube experiment, resulting from the given neutrino 
spectra; the new curves all fall less steeply with energy than before. Above 1 TeV, the new curves are higher by up 
to an order of magnitude. Below 1 TeV, the new curves are slightly lower for "All Muons" (meaning a trigger on all 
muon tracks), and an order of magnitude lower for "Muons (Contained Vertex)". Thus, the statistical uncertainties 
will be improved when muons above 1 TeV are used, which is preferred, since this is well matched to the 1 TeV 
threshold required for resolving showering events. (Previously, we had argued that using a muon threshold of 100 
GeV was statistically beneficial, at the cost of some model dependence; this argument is no longer necessary.) As 
before, contained vertex muons offer a more robust neutrino energy estimate than do the throughgoing muons that 
dominate the all-muons sample, although the statistics of the all-muons sample is better. The rate of the contained 
events is still about 1/10 of the all-muons sample, but now the absolute statistics of both are more favorable. 

Figures 7 and 8 reveal graphically the most important results of our paper. They show the expected ratios of muon 
events to showering events (all events other than muon-track events) , as a function of the v^, fraction of the neutrino 
flux. The IceCube experimental efficiencies are folded into these figures. Thus, these figures provide the map from 
the incident flavor fraction onto a direct experimental observable, and vice versa. (We assume v-^ interchange 

symmetry, which is supported by inferences from neutrino oscillation results.) The qualitative message of these figures 
is the same as before: an IceCube measurement of the muon/shower ratio can determine the flavor content of the 
incident neutrino beam with sufficient precision to study neutrino physics and astrophysical dynamics. 

While the normalizations of the new curves are quite different, the statistics are increased, so that the final conclu- 
sions are quite similar. For example, consider an incident fraction of 1/3, as expected from pion decay, and note 
the predicted muon/shower ratio. Using the dashed and dotted lines, corresponding to the uncertainties, note also 
the range of different incident v,, fractions that would give the same muon/shower ratios. In this standard case, the 
Vf. fraction would be deduced to be in the ranges 0.25-0.40 and 0.14-0.50 for the two fluxes in the new Fig. 7 (100 
GeV muon threshold), and 0.24-0.41 and 0.12-0.52 for the two fluxes of the new Fig. 8 (1 TeV muon threshold). 
The ranges quoted in the published paper were similar: 0.26-0.37 and 0.18-0.44 for the lower muon threshold, and 
0.22-0.42 and 0.09-0.61 for the higher muon threshold. The muon/shower ratios in the flnal column of Table I should 
be changed from the original values of (9, 1.5, 40, 14) for z^e-fraction (0.33, 0.75, 0, 0.2), to (3, 1.0, 5, 3.7) for the same 
i/g-fractions. Thus, while the normalizations of the curves have changed, the precision with which the flavor ratios 
can be measured has not. 

In summary, with the corrections to our paper, the measurements of neutrino spectra and flavor ratios are made 
easier, especially above 1 TeV. The other results of the paper, including the detection probabilities for different 
neutrino flavors, are unchanged. 

We are grateful to Shin'ichiro Ando for pointing out that our results were in error, and for subsequent discussions. 
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FIG. 9: The replacements for Figs. 1 and 2. 
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